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Abstract 

We obtain asymptotic expansions for Toeplitz determinants corresponding to 
a family of symbols depending on a parameter t. For t positive, the symbols are 
regular so that the determinants obey Szego's strong limit theorem, li t — 0, 
the symbol possesses a Fisher-Hartwig singularity. Letting t — >■ we analyze 
the emergence of a Fisher-Hartwig singularity and a transition between the two 
different types of asymptotic behavior for Toeplitz determinants. This transition 
is described by a special Painleve V transcendent. A particular case of our result 
complements the classical description of Wu, McCoy, Tracy, and Barouch of the 
behavior of a 2-spin correlation function for a large distance between spins in the 
two-dimensional Ising model as the phase transition occurs. 

1 Introduction 

Consider the Toeplitz determinant with symbol f{z) G L^(C), where C is the unit 
circle: 

1 /■2n 

= det(/,_fc)^;io, /, = -j^ f{e'')e-^^'d9. (1.1) 

We are interested in the behavior of Dn as n — )• oo. 

If In f(z) is sufficiently smooth on the unit circle (in particular, f{z) is never zero 
for z G C and has no winding around the origin) so that ln/(z) G L^{C) and the sum 

J] |fc||(ln/)fc|2, (ln/), = — y lnf{e'')e-'>''de, (1.2) 

k=—co ^ 

converges, then the asymptotics of Dn are given by the strong Szego limit theorem 
[221 [m [25]: 

lnDn = — ln/(e*^)de + J^A:(ln/)fc(ln/)_fc + o(l), as n ^ oo. (1.3) 

^ -^^ k=i 

However, one often encounters a situation where the symbol possesses so-called 
Fisher-Hartwig singularities. In the case of only one such singularity, located at z = 1, 
the symbol has the form: 

f{z) = Iz-lp^z^^e-^'^e^^") = (2-2cos6')°e^^(''-")e^("'''), for < ^ < 2tt, (1.4) 
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where V{z) is a sufficiently smooth function (see pH]) on the unit circle. The singularity 
at z = 1 combines a jump-type (for a = 0, /3 / 0) and a root-type singularity (for /3 = 0, 
Q 7^ 0). For this symbol the sum ()1.2p diverges, and therefore Szego's theorem does not 
hold. The asymptotics for the Toeplitz determinant are given instead by the expression 

oo oo oo 

In = nVo + Y, ^^^V-k - (a - /3) ^ 14 - (a + /?) V.u 

k=l k=l k=l 

,r2 a2^^ G(l + O + /3)G(1 + a - /3) , , ,^ . 

+ (a^ -/3^)lnn + ln^ ^^^'^ — + o{l), as n ^ oo, (1.5) 



if 



a±/?/-l,-2,..., (1.6) 



with 
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27r 







Here G is Barnes' G-function, which is an entire function having the properties: G{z + 
1) = r(z)G(z), where T{z) is Euler's T-function, and G(l) = 1, G{-k) = for A; = 
0,1,2,.... Note that if V{z) = 0, there exists an explicit expression for D„ with symbol 
\z — l\'^"z^e~^'^^ in terms of G- functions [9l [TO]. 

Suppose now that a symbol depends on a parameter t {f{z) = f{z; t)) so that when 
t > the symbol is "regular", i.e. Szego's theorem holds for Dn{t), while at t = the 
symbol has the form (II. 4p . The purpose of the present paper is to study the transition 
from (jl.3p to (|1.5p as t — )■ 0. Namely, consider the following symbol 



f{z) = {z- e^Y^\z - e-*)°-/32-°+/3e-*^("+/3)e^(^), q ± /3 / -1, -2, (1.7) 

where t > is sufficiently small and a, /3 G C with Re a > — ^. We further assume that 
V{^z) is analytic in an annulus containing the unit circle and write it there in terms of 
its Fourier series 

+ 00 

V{z)= VkzK (1.8) 

fc=— oo 

We define the powers in (jl.7p with arguments between and 27r. With this choice 
of branch cuts, /, and moreover In/, is analytic in C \ ([0,e~*] U [e*,+oo)) and, in 
particular, on the unit circle for i > 0. Therefore, for any fixed f > 0, the asymptotics 
of the Toeplitz determinant -D„(t) are given by ()1.3p . Calculating the Fourier coefficients 
(ln/)fc, we obtain 



In £)„(*) = nt{a + 13)+ nVo 



+Y^k 

k=l 



Vk-{a + (3) 



g tk-\ 



V-k -{a -13) 



e 



+ o(l), as n — co, i > 0. 

(1.9) 



For t = 0, the symbol reduces to (jl.4p (with analytic V). Therefore, for t = 0, the 
asymptotics of Dn{t) are given by ()1.5p . In the present paper we describe the transition 
from (11.91) to (11.51) when t decreases to 0. 



2 



The paper also sets the stage for analysis of various other transition asymptotics 
for Toeplitz determinants, such as 2 singularities approaching each other; emergence of 
an arc of the unit circle where the symbol / = from 2 jump-type singularities at the 
ends of the arc, etc. 

Our analysis explains to some extent the question of connection between Painleve 
tau-functions and Toeplitz determinants which was noticed before: see [7j for a discus- 
sion and references. Historically and as the most prominent example, this connection 
appeared in the study of 2-spin correlation functions in the 2-dimensional Ising model. 
We discuss this in the following section. 



1.1 Application: the two-dimensional Ising model 

Transitions between Szego weights and Fisher-Hartwig weights arise for example in 
the theory of solvable two-dimensional statistical models and one-dimensional Heisen- 
berg spin chains. Recall the two-dimensional Ising model solved by Onsager (see, e.g., 
|32^ [29j ) . In this model a 2Ai x 2J\f rectangular lattice is considered with an associated 
spin variable ajf^ taking the values 1 and —1 at each vertex {j, k), —M. < j < — 1, 
—N < k < N — 1. There are 2^^^ possible spin configurations {a} of the lattice 
(a configuration corresponds to values of all 0"^^ fixed). We associate with each con- 
figuration the energy of the nearest-neighbor coupling (imposing the cyclic boundary 
conditions on the lattice) 

M-i N-i 

E{{<y}) = - ^ ^ {-fiajkCTj k+i + l2<y jk(T j+ik) ■, 71,72 >0. (1.10) 
j=-M k=-N 

The partition function at a temperature T > is equal to 

Z(T) = ^e-^«->)/T^ (1,11) 
{^} 

where the sum is over all configurations. A remarkable feature of this model is the 
presence of a thermodynamic phase transition in the limit of the infinite lattice at a 
certain temperature Tc depending on 71, 72. 

Define a 2-spin correlation function by the expression 

<fT00fTn„ >= ^ lim -L- VcT00CTnne-^({"»/^. (1.12) 

For large n, this function measures the long-range order in the lattice at a tempera- 
ture T, which determines magnetization. Indeed one can show that the spontaneous 
magnetization M is given by the expression 



M = / lim < aooann >■ (1-13) 

It is a remarkable fact that the 2-spin correlation function is a Toeplitz determinant 

< aooann >= e^'/^Dn{t), f{z;t) = {z- e')-'/\z - e-')'/^z-'/^e'-/^ (1.14) 

where 

e* = sinh sinh , (1-15) 
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and the branches of the roots are chosen with the arguments from to 2tt. The symbol 
in ()1.14p has the form ()1.7p with 



a = 0, /3 = -2' V{z) = 0. (1.16) 

The critical temperature Tc is defined by the condition that t = 0. 

For T < Tc we have t > 0. Therefore the strong Szego limit theorem (jl.9|) holds, 
and using the elementary identity 

oo 

e-^^'/k = - ln(l - e-2*), (1.17) 

k=l 

we rederive the well-known result 



e«*/2z)„(t) = (i_e-2*)V4(i + o(i)) 

1/4 

271 272 ' 
smh — — smh — — 
T T 



(1 + 0(1)), as n 00, T < Tc. (1.18) 



The correlations tend to a constant as n — )• 00: the model exhibits the long-range order. 
Notice that by ()1.13p . the asymptotics (jl.lSp imply the existence of the spontaneous 
magnetization for T < Tc and the famous power-law decay of the magnetization M ~ 
Const • {Tc - Tf/^ as T / Tc. 

At T = Tc, we have t = 0, and therefore a Fisher-Hartwig singularity with the 
parameters a = 0, /3 = —1/2 appears sX z = 1. In this case, (jl.Sp holds and we obtain 
(cf. |32l|29]) 

DniO) = ^ \' ' (1 + 0(1)), (1.19) 

so the correlations decrease as n~^^^, n — )• 00: the long-range order is destroyed. Note 
that as V{z) = 0, there is an explicit expression for -Dn(O) (cf. the remark following 

2 \ n n—l / ^ \ k—n 



o»(o) = ^) n 

^ ^ k=\ 



For T > Tc we have t < 0. The symbol still has a singularity at z = 1 but now 
with the parameters a = 0, /3 = — 1. This is the situation of a degenerate type of a 
Fisher-Hartwig singularity, and (jl.Sp does not hold in this case as one of the G-functions 
vanishes. Calculations [32] show an n~^/^e"* decay of the correlations (jl.l2p as n — )• c«. 
There is no long-range order and M = 0. 

The transition T — )• for large n was studied in |36l [30l [33] (for a more general 
correlation function < aooanm >)■ The authors took n — )■ 00 with x = n(e^* — 1) fixed 
and found, in particular, that in this limit 

n^/^ < aooann >^ F{x), (1.20) 

where F{x) is given in terms of a solution to Painleve III equation (reducible to Painleve 
V: see equation ()1.28p below). Moreover, in [30], McCoy, Tracy and Wu evaluated 
the connection formulae for this Painleve III function and showed that the limiting 
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behavior of F{x) as x — >■ oo and as x — )• formally matches the asymptotics (jl.l8p 
and (jl.l9p . respectively. The matching with (|1.19p was, however, checked only up 
to a multiplicative constant. A more detailed evaluation of the small x behavior of 
the function F{x) was carried out later by Tracy in [M], and reproduced exactly the 
constant in the critical point asymptotics (jl.lOp . The calculations of [30] and [M] 
were based on an alternative representation of the function F(x) as an infinite series of 
integrals and were rather involved. 

The results of the present paper, namely a particular case of Theorem 11.41 below, 
fully describe (up to a possibility of a finite number of poles: see below) the transition of 
the correlation function ()1.12p for large n from T < Tc to T = Tc, that is the transition 
from (jl.lSp to (jl.lOp . Note that we do not rely on (jl.20p as the parameter x in our 
analysis is not necessarily fixed, in fact, our asymptotics for the correlation function are 
uniform in the whole range x £ [0, oo) (away from a finite number of positive points). 
More precisely, our asymptotics as n — >• oo are uniform for all T G [Ti, Tc] (away from 
a finite number of positive x's) for some Ti < Tc with Tc — Ti sufficiently small. 

The description of the Ising double scaling theory which we obtain as a particular 
case of Theorem 11.41 is in agreement with the classical results of |36[ [30l [M] (see also 
Remark 1 1 . 71 below) . Furthermore, by obtaining the uniform asymptotics for the whole 
transition range of temperatures T < Tc, we complement the analysis of this case of 
the Wu-McCoy-Tracy-Barouch scaling theory. 

Another example where Theorems 11.11 11.41 below can be applied is the so-called 
emptiness formation probability in a Heisenberg spin chain [20j . 

1.2 Statement of results 

Consider the second order ODE 

cP'a\^ ( do ^ f da\'^ „ d(j\ 



dx"^ J \^ dx~^ \dxj ^ dx j 
This is the Jimbo-Miwa-Okamoto cj-form |25l |26] of the fifth Painleve equation 



1 1 \ ^ I (u-l)2/ B\ Cu u{u + l) 

— + ]ul--u, + ^ ^ {Au + -]+ — + , 

2u u — 1 J X V u J X u — 1 



with the parameters A,B,C,D given by 

^ = l(a-/3)2, B = -l{a + /3f, C = 1 + 2/3, ^ = - J- (1-23) 



In the following theorem, we give the asymptotic expansion for the Toeplitz deter- 
minant with symbol ()1.7p as n — )• oo which is valid uniformly for < t < to- Our 
asymptotic expansion interpolates between Szego and Fisher-Hartwig asymptotics. 

Theorem 1.1 Let a G M, a > — ^, /? G iM. Let f be defined by |j. ?[ ) and consider the 
Toeplitz determinant Dn{t) defined by il.l]) corresponding to this symbol. The following 
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asymptotic expansion holds as n —)• oo with the error term o(l) uniform for < t < to 
where to is sufficiently small: 



InDnit) = uVq + (a + /3)nt + 



fc=i 



Vfc - (a + /3) 



e 



- (a - /3) 



^^""^ G(l + 2a) '^ + n{2nt) + o{l), (1.24) 

where G{z) is Barnes' G-function, and 

ni2nt)= ^ dx + {a'^ - /3^) In 2nt. (1.25) 



T/ie function cr{x) is a particular solution to the equation ^.21\) which is real analytic 
on (0, +oo), and has the following asymptotics for x > 0: 



a{x) 



^2 _ ^2 ^ 9^{x - x^+^'^Cia, /3)}(1 + 0(x)), X ^ 0, 2a ^ Z 
a2 -/32 + 0(x) + 0(x^+2") + 0(xi+2"ina;), x ^ 0, 2a G Z (1.26) 



l+2a„-x -1 



X ^ e 



r(a-/3)r(a+/3) 



;i + 0(i)), x^+oo, 



with 



^ r(l + a + /3)r(l + a-/3) r(l - 2a) 1 
L^l«,PJ r(l-a + /3)r(l-a-/3)r(l + 2a)2l + 2a' ^ ^ ^ 

where T{z) is Euler's T-function. 

Remark 1.2 Later on, we will construct (t(x) explicitly in terms of a Riemann-Hilbert 
problem. 

Remark 1.3 With increasing effort, one can calculate more terms in the expansion 
(I1.24P using our approach. 

The function a = cr{x; a, /?) is defined for x E C with a cut from zero to infinity. It 
is analytic in the cut plane apart from possible poles. Asymptotics (|1.26|) imply that 
there are no poles for x positive and sufficiently large. Hence the number of possible 
poles of ct(x) on (0,+cx)) is finite. We show below that for a > — ^ real, (3 imaginary, 
there are no poles on the real half-axis (0, +00). Therefore we took the intervals of the 
real line as a path of integration in (ll.25p . For (3 arbitrary. Re a > —1/2, a similar 
result holds, but we have to choose a path of integration in the complex plane avoiding 
possible poles which we denote {xi, . . . , xe}- Namely, we have 

Theorem 1.4 Let a, /3 G C with Re a > — ^, a it /? 7^ —1, —2, . . . , and let ss denote 
a sector —tt/2 + 5 < argx < 7r/2 — 5, < 5 < -k/2. Let f be defined by |j. Tj j and 
consider the Toeplitz determinants Dn{t) defined by corresponding to this symbol. 
There exists a finite set {xi, . . . ,X£} E ss (with I = £{a,(3,6) and Xj = Xj{a,(3) ^ 0) 
such that the expansion ^1-24^ holds uniformly for t G ss,\t\ < to (with to sufficiently 
small) as long as 2nt remains bounded away from the set {xi, . . . ,xe}. The function 
^} is defined by il.25\) . where the path of integration is chosen in ss, connecting with 
2nt and not containing any of the points {xi, . . . ,X(}. Moreover o"(x) solves the ODE 
hl.21\) and has the asymptotics in the mentioned sector given by [T7i 
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Remark 1.5 It follows from the representation (jl.24p that the residue of ^cr{x) at 
each of its poles in the sector — 7r/2 < argx < 7r/2 is an entire number. Different 
choices of the integration contour in ()1.25p correspond, in general, to different branches 
of In Dn{t). If a and /3 are such that has a pole x, the determinant Dn{t) is zero 

at 2nt = X (up to an 0{l/n) error term). 

Remark 1.6 From now on, we will always consider t > for simplicity. The extension 
to t G S5 is straightforward. 



Remark 1.7 In the example of the Ising model discussed in the previous section, we 
have « = 0. B = —1/2 (see I1.16"]) . and equation (jl.2ip becomes 




(1.28) 



This is exactly the equation which was obtained in jt26j for the function 



ax) = x^\nF{x)--^, (1.29) 

where F{x) is the right hand side of the Ising double scaling limit (jl.20p (see also 
equation (4.16), with r = x/2, in [29]). It follows immediately from our main result 
(I1.24P that the function C,{x) in |26j and our function (j{x) coincide: 

C{x) = (t{x). (1.30) 

Thus the application of our Theorem 11.41 to the Toeplitz determinant (I1.14p yields the 
complete analysis (up to the question of existence of a finite number of positive poles 
Xj) for T < Tc of the Jimbo-Miwa Painleve V version of the Wu-McCoy-Tracy-Barouch 
scaling theory for the 2D Ising model. 

From the expansion (11.24p . we can recover the Fisher-Hartwig asymptotics for 
lnL>„(0). Let t 0, and n fixed in (fL2l) . Then, using (fr25|) and (fl^ . we ob- 
tain that n{2nt) = (q^ - p"^) ln(2nt) + o(l) if Re a > -\. Substituting this into (fOi]) 
and recalling (jl.lTp . we obtain (jl.Sp . 

The expansion ()1.24p should also be consistent with the Szego asymptotics for t 
fixed. We see immediately that the 0{n) term gives, for a fixed t, the corresponding 
term in the Szego asymptotics. Consistency of the 0(1) terms, however, yields an 
interesting identity involving the Painleve function (j{x) via (I1.25P : 

o^^ ^ 1 G'(l + a + /3)G(l + a-/3) ,^ 
1.3 The Painleve V Riemann-Hilbert problem 

We can say more about the function (t{x) than we did in Theorem ll.lt we can construct 
it explicitly in terms of a Riemann-Hilbert (RH) problem. Consider the contour F = 
yj^j^iTj in the complex plane (see Figured]), with 

r4 = ^ + e^^M+, F5 = (l,+oo), F6 = (0,l), 
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with Fi, . . . oriented towards infinity and Fg oriented to the right. By a standard 
convention, the "+" side of the curve is on the left as one faces the direction of the 
curve's orientation. 

Let Rea > — ^ and consider the following RH problem for * = '^{C; x, a, P). 

RH problem for * 

(a) * : C \ F ^ is analytic. 

(b) ^ has continuous boundary values on F \ {O, ^, 1} , and they are related as follows, 

*+(C) = *-(C)(^Q 1 j' forCGFi, (1.32) 

*+(C) = *-(C)(_^J(a-^) °), forCGF^, (1.33) 

^+(C) = ^-(C)(^.,(1-^) J), forcers, (1.34) 

*+(C) = ^-(C) ' -, , forCGF4, (1.35) 



vI/+(C) = ^-(C)e'"'^"=\ forCGTs, (1.36) 

*+(C) = *_(C)e-"("-^)'^3^ fo, ^ g Te, (1.37) 

with £73 = 

(c) * has the following behavior as C ^ oo (for some matrices Ci = Ci{x,a, /3), 
C2 = C2{x,a,P)), 



(do) As C 0, 



*(C) = o(|':|^ l^l-H^I- (1-39) 
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(dl) As C ^ 1, 



^(c) = of!^ ^ll,]- (1-40) 



IC-il 



Furthermore ^' is bounded near ^. 
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The RH conditions imply (by a standard argument) that the determinant of the solution 
^ (which is, if it exists, unique) is identically equal to 1, and consequently we have 
using (11.381) that tr Ci = 0. Let us denote the matrix elements of Ci by 



C^(x) 



q{x) r{x) 
t{x) —q{x) 



Define the functions v and u in terms of the matrix elements of Ci: 

^(^) = — ~ ~ xr{x)t{x), (1-41) 

xt 

u{x) = l + --: -— T--. (1.42) 

^ ^ {2l3 + l-x)t{x)+xt'{x) ^ ^ 



We will show in Section [4.31 below that 

r+oo 

a{x) = / v{C)dC (1.43) 



is the function appearing in Theorem 11.11 and Theorem 11.41 The RH problem for ^ 
is a special case of the RH problem associated to the fifth Painleve equation, see e.g. 
MM- 

We prove the following. 

Theorem 1.8 Let a,/3 E C and Re a > — ^. 

(i) The RH problem for ^ is uniquely solvable for all x > except possibly for a finite 
number of positive x-values. We denote the x-values for which the RH problem 
is not solvable by {xi, . . . with /3) and k = k{a, f3). 

(a) Iflma = and Re/3 = the RH problem is solvable for all positive x-values. 

(Hi) The function v defined by |i..^i[ ) is analytic in (0, +oo) \ {xi, . . . , Xk}, and solves, 
together with u defined by the system 

xux = xu- 2v{u - if + {u- l)[(a - /3)n - /? - a], (1.44) 

V 

xvx = uv[v — a + f3] (f — 13 — a). (1-45) 

u 

(iv) The function v has the asymptotics given by 



v{x) 



-^^2^{1 - (2a + l)x2°C7(a, /?)}(! + 0{x)), x ^ 0, 2a ^ Z, 

0(1) + C'(x2") + C)(3;2"lnx), x^O, 2a e Z (1.46) 



X ^'""e ■^w:r-^3w:-rm (1 + O (^)) , x +oo. 



-l+2a„-x -1 

r(a-/3)r(a+/3) 



where C{a,/3) is defined in fL2^. 
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In addition, we have 
r+oo 

/ v{x)dx = a2 _ ^2 
Jo 

i/Ima = and Re /3 = 0. In the general case, a,/3 G C, Rea > equation 
(Ljl^ holds up to addition of 2iiim, m G Z, with the path of integration avoiding 
{xi, . . .,Xk}. 

Part (iii) of the theorem follows from a standard Lax pair argument and was proved 
in [18', TB] for a slightly different but equivalent RH problem. That proof applies to 
our RH problem as well, and implies moreover that the RH solution is meromorphic 
in X for x G C \ {0}. We will come back to this in Section HI We prove part (iv) 
by performing the Deift-Zhou steepest descent analysis for the RH problem. This 
asymptotic analysis also implies the solvability of the RH problem for large x and 
small x^ and by meromorphicity in x, this leads to the statement (i). We prove part 
(ii) by applying the technique of a vanishing lemma to the RH problem for ^. 

Remark 1.9 The system (ll.44p - (11.45p is related to the Painleve V equation: eliminat- 
ing f , we easily verify that u solves the Painleve V equation p.22p - ()1.23p . Asymptotic 
expansions as x — )• and as x — )• oo for various solutions to the fifth Painleve equation 
and the system (jl.44p - (jl.45p were obtained in several works, see e.g. [T| l2t l3l [25 } [3 H [33] . 
The solution v which is of interest to us decays exponentially at +oo, is integrable near 
if Rea > — |, and it has no poles on (0, +oo) if a > — ^ G M and (3 G iM. Note 
that the asymptotics of a{x) (f06]l fohow from (fL46]l and (fLiTP by (fL43]l . Although 
we obtained the connection formulae ()1.26p as a by-product of our analysis, we would 
expect that these asymptotics can be found in the general list of connection formulae 
for the fifth Painleve equation obtained in [H EJ [3] . 

Outline of the paper 

The proofs of Theorem 1 1 . 1 1 and Theorem 11.41 are based on a well-known connection be- 
tween Toeplitz determinants and orthogonal polynomials on the unit circle. In Section 
O we obtain a differential identity for lnZ)„(t) in terms of the polynomials orthogonal 
on the unit circle with weight f{z). In Section [3l we obtain large n asymptotics for 
these orthogonal polynomials from a RH problem. The asymptotics will be given in 
terms of a model RH problem which we study in detail in Section HI where we also 
give a proof of Theorem 11.81 In Section [5l we use the previously obtained asymptotics 
for the orthogonal polynomials and the results of Section S] to integrate the differential 
identity for lnD„(t), which leads to Theorem 1 1 . 1 1 and Theorem 11.41 

Throughout the paper, we choose the branches of logarithms and roots correspond- 
ing to arguments between and 27r, unless stated otherwise. 

2 RH problem for orthogonal polynomials and a differen- 
tial identity for the Toeplitz determinants 

Our analysis is based on a classical connection between Toeplitz determinants and 
orthogonal polynomials. Assume that for some n > Dn,Dn+i ^ 0, and define a 
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polynomial (pni^) in terms of the Fourier coefficients of f{z) as follows: 



4>n{z) = 



n+1 



fo /-I 

/l fo 

fn-1 /n-2 
1 Z 



/-n+1 
/-I 



(2.1) 



The leading coefficient of ^„ is then equal to 



y Dn+i 

There holds the orthogonality relation 



j = 0,l,...n, 



(2.2) 



(2.3) 



where C is the unit circle oriented in the counterclockwise direction. Similarly, let 
(pniz) be defined by 



1 



n+1 



fo f-l 
fl fo 



/-n+1 1 
f-n+2 Z 



(2.4) 



fn fn-l ■ ■ ■ fl Z 

Then (pn has the same leading coefficient x„ as and 

[ (t)n{z'^)z^ fiz)^ = Xn^Sjn, j = 0,l,...,n. 

l-K Jc IZ 



(2.5) 



If Dn 7^ for n = 1, . . . (and we set Do = 1, (l)o{z) = 0o(^) = ^/VDi), the sys- 
tem of polynomials ^„ and <^„, n = 0, 1, . . . exists and can be characterized by the 
orthonormality relations 



TT- / (t>k{z)4>m{z ^)f{z)— = 5km, 

27r Jc tz 



k,m = 0,1, ... . 



(2.6) 



If the symbol / is positive on the unit circle C, it is a classical fact (which follows, 
e.g., from the representation of a Toeplitz determinant as a multiple integral) that 
Dn{f) > for all n > 0, and the system of orthogonal polynomials exists. 

Assume that D^-i, D^, 7^ 0, t > 0, and define the function Y{z; n) as follows 



Y{z) = 



-1 r 0n(g) /m 



Xn^(t>n{z) Xn^ Jc 

-Xn-lZ--'k-l{z-') -Xn-l/c%S^^ 



(2.7) 



Then Y is the unique solution of the following RH problem with a jump on the coun- 
terclockwise oriented unit circle C. 
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RH problem for Y 

(a) y : C \ C ^ C2>^2 ig analytic. 

(b) Y+{z) = Y^{z) (J ^^Y^^^) ' for ^ e C. 

(c) Y{z) = {I + 0{l/z)) (^J , as z ^ oo. 

A general fact that orthogonal polynomials can be so represented as a solution 
of a RH problem was noticed in [T7| (for polynomials on the line) and extended for 
polynomials on the circle in [3]. 

In the next section we will show that the RH problem for Y{z; n, t) is solvable (and 
therefore the orthogonal polynomials exist and the coefficients Xn are nonzero) for all 
n larger than some no{a,/3) provided 2nt is bounded away from a certain finite set 
of points (in particular, see Proposition 12.11 below. 7^ 0). The number no(a,/3) is 
bounded for a and /3 in a bounded set. 

Our next aim is to express ^ In Dn{t) in terms of the entries of the RH solution Y. 
We prove the following. 

Proposition 2.1 Let t > and n G N. Suppose that the RH problem for Y{z;n,t) is 
solvable. Then Dn / 0, and the following differential identity holds: 

± lnDr.{t) = -{a + /3)e'(Y-^^) (e*) + (a - /3)e-* f y-i^) (e"*). (2.8) 
dt \ dz J 22 V dz J 22 

Proof. We will follow the approach of [24J. Let us start with the expression 

Dn{f) = <lei{I-Kn), (2.9) 
where Kn is an integral operator acting on L^{C) with kernel 

, (z/z')'' -11- f(z') 
K4z,z')=^-^^^ (2.10) 

Z — Z' ITTI 

This fact is easy to verify by considering the matrix expression for Kn in the basis {z'^}, 
k = —00, . . . , c«. 

As follows from the theory of "integrable" Predholm operators (see, e.g. [24J), the 
solvability of the RH problem for Y{z; n, t) implies that the operator 1—Kn is invertible. 
Therefore, £»„(/) = det(/ - K^) / 0. 

Consider 

I InDnif) = |tr ln(/ - Kn) = -tr (/ - Kn)~^'^. (2.11) 



Since 

^^-^ " + ^e* + ^^e-M/, (2.12) 



dt 

we have 
dK 



kW-k(,'), (2.13) 
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where 



and 



% ^ ^ Z Z/T^% 

Here 

A.(...-)^i^^('/^')°7V (2.16) 

^ ' ^ z'-e*27ri z-z' ^ ^ 

The reason to single out Ai will soon become clear. By a residue calculation, we obtain 



. — + " ^ ' ^ 



27rf z — e* 
We can now rewrite (j2.15p as follows: 

z') = (Ai(/ - K^)){z, z') + {a + p)e'l^l^^^Mpl^^, (2.17) 

Im [z — e^)[z' — e^) 

Defining the following 2-component vectors 

we can write (j2.17p in the final form: 

K^^\z,z') = {K^{I-Kn)){z,z') + {a + P)e\e-^'h{z)Uz')-h{z)^^^^^ (2.18) 
On the other hand, let us define Fj = (/ — Kn)~^ fj- Then 
F,{z) f Kn{z,z')^..^, 



z — & Jc z — e 
Noting that 

J, hzfg{z') 



Kn{z,Z 



we can write the above equation in the form 

^'^'^ /■ K^{z,z')^^dz' + / J{zfgiz')EMldz' = f,{z). 



z — e* Jc ' z' — e'- Jc ' z' — e 

Applying (/ — Kn)^^ to both sides, we obtain 

2 

z — e 



2 ^ 1 

Vm,fc(e*)Ffc(z) = -Fj{z), j = 1,2, (2.19) 

^ — ' z — e'- 

k=l 

where F^ = {I - K^)'^ fk and 

mjkie') = 5jk- / Fj{z)gk{z) -, j,A; = l,2. (2.20) 

Jc z — e 
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Thus we have for the 2-componeiit vector (detm = 1, see 



^(^) = :— I"^ (e)^(^) = — — 7 Z7_ , _..i7J • (2-21) 



As is shown in ([23], Eq. (2.16) up to a different notation), the matrix m is related to 
y. For \z\ > 1, 

y— m2lZ — 77122 m22Z J 

Using the definition F = {I - Kn)~^f, equations (IXTH|) . (IX^ . (IX^ . and the fact 
that 

e* f ndz' 



trAi = (a + /3)— / — ^ =-n{a + P), 

we easily obtain 

tr ((I - K„)-^k(^)) = (a + /3)e*(yn(e*)i;'22(e*) " Y2i{e')Yl,,{e')), (2.23) 

where ^.'(e*) stands for the derivative of Y(yz) w.r.t. z evaluated at z = e*. 

(2) (2) 
Let us now compute the contribution of Kn ■ First, write Kn in the form 



^(2) . a-/3 (z/zQ" - 1 1 - /(zQ a-/3 (^/z^)" - l l ^ 
' = A2 : 7 : — : e , = — 7 : — rC 



We then obtain as above for Kn ^ that 

Ki^\z, z') = (Mil - K„))(z, z') + (a - (2-24) 
Defining the new vectors 

fiz) = J^, giz) = J^, F, = {I-K)~'fk, 

z — e ^ z — e ^ 

we can write (I2.24p in the form: 

k(2)(2, = (A2(/ - Knmz, z') + {a- (3)e-\e-^'f2{z)gi{z') - fi{z)gi{z')) (2.25) 
and obtain as above 

p. X _ 1 ( "i22(e-*)Fi(z) - mi2(e-*)F2(z) \ 

^ ^ z- e-* V-"^2i(e-*)Fi(z) + mii(e-*)F2(z) J ' ^^'^"^ 
For |z| < 1, the matrix m is related to Y by the formula preceeding (3.21) in [24] : 

y(z)=f "^""" + "^12 "^^0. (2.27) 
\-m2lZ"- - 17122 -m.21/ 

Using the definition F = {I - Kn)'^!, equations ([225]), ([226]), dZ^Z]), and the fact 
that tr A2 = 0, we finally obtain 

tr ((/ - Kr^r'KP) = {a- P)e-\Yn{e-')Y^,,{e-') - Y2i{e-')Y^,,{e-')). (2.28) 

Expressions (j2.23p and (|2.28p imply by (j2.13p the statement of the Proposition. □ 
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3 Asymptotic analysis of the RH problem for orthogonal 
polynomials 

In this section, we apply the steepest descent method of Deift and Zhou [13] to the RH 
problem for Y. We follow the strategy used in [3| for orthogonal polynomials on the unit 
circle with a special weight function, and used in [12\ TlSj for orthogonal polynomials on 
the real line with respect to a more general weight function. The most important new 
feature here is the construction of a local parametrix near 1 which involves a Painleve V 
RH problem. We will obtain asymptotics for y as n — t- cxd, and therefore, asymptotics 
of the r.h.s. of (12.811 in terms of Painleve V functions. 



3.1 Normalization of the RH problem 

Define 

n^) = {::r:~ ' (3-1) 

as \z\ < 1, 




with Y given by (12. 7p . Then T satisfies a RH problem normalized at infinity. 



RH problem for T 

(a) T ■.C\C^ C2X2 is analytic. 

(b) T+(z) = T.{z) , for z G C. 

(c) T{z) = 1 + 0{l/z) as 2 ^ oo. 

The diagonal elements of the jump matrix for T oscillate rapidly on the unit circle if n 
is large. The next transformation turns the oscillatory behavior into exponential decay 
on a deformed contour. 

3.2 Opening of the lens 

Note that one can factorize the jump matrix for T as follows. 



Mz) 



f{z) 



1 0\ / f{z 

ij \-fiz)-^ ; \z^f{z)-^ 1 



(3.2) 



To make use of this factorization, consider the three counterclockwise oriented closed 
curves as shown in Figure [2j Let us write 



J^(') = [,-nf(^,yi ij> (3-3) 

M^)=(.n,L-l ?V (3.5) 
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and define 



S{z) 



T{z), in regions ili and 

T(z)Ji, in region 0,2, (3-6) 
^T(z)J^^, in region fl^, 



with 0,1,0,2,^3 as indicated in Figure [2j Note that / is an analytic function in C \ 
([0,e~*] U [e*,+cxD)), and we define Si, S2 so that the branch cuts are located in the 
regions Oi and O4. As we can have i — t- 0, it is inevitable that Oi and O4 approach 1. 
We choose Si and S2 to go through 1 as in Figure [21 The function S{z) satisfies the 
following RH problem. 

RH problem for S 

(a) S : C \ (Si U C U S2) ^ C^xs is analytic. 

(b) S+{z) = S-{z)Jk{z), for z G Sfc, A: = 1, 2, 
Sj^{z) = S^{z)Jn{z), for z £ C. 

{c) S{z) = I + 0{l/z), asz^oo. 

Note that the jump matrices Ji and J2 tend to the identity matrix on their respective 
contours Si and S2 as n — )• 00 except near 1. 

We need to construct a parametrix dealing with the jump condition on the unit 
circle and a local parametrix near 1. 

3.3 Global parametrix away from 1 

Ignoring the jumps on Si and S2 and a neighborhood of 1, we consider the following 
model problem. 

RH problem for A'^ 

(a) iV : C \ C ^ C2x2 is analytic. 

(b) N+(z) = N^{z)Jn{z), for z € C. 
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(c) N{z) =1 + 0{l/z), as z oo. 
This problem is easily solved explicitly: 



N{z) 



D{zY^ ( ° I , for |z| < 1, , , 

^ ' \-l Qj' ' ' (3.7) 

^D{zYi, for \z\ > 1, 



where the (Szego) function D(z) is analytic and nonzero in C\C, tends to 1 as z ^ oo, 
and satisfies the jump condition D^{z) = D^{z)f{z) for z & C. It is easy to verify that 



_ ;(^-e*)"+/^e— (-+/^)exp(Er=o^fc^') > fo^ 1^1 < 1, 

^ {z - e-'r^^^z-P exp (- Efci-oo ykz"") , for \z\ > 1. ^"^'^^ 



3.4 Local parametrix near 1 

For < i < io with to fixed but sufficiently small, we will now construct a parametrix 
P satisfying the same jump conditions as in a neighborhood [/ of 1 of a sufficiently 
small fixed radius and a matching condition with on the boundary dU . 
Assume that ^(C) solves the RH problem of Section [1.3| and define 

$(A;x) = e-i^'^'x-f^^'-'-^i- + -; x)G(A; x)^'"='e^^("-'^)'^3^ for iblmA > (3.9) 

X 2 

respectively, with 

G{\;x) = (A + |)"("~^)(A - |)"+/3e^e-"^("-^), x > 0, (3.10) 

where G is analytic in C \ ((— oo,— |] U [|,+oo)). We choose — vr < arg(A + |) < vr 
and < arg(A — |) < 2it. It is straightforward to check that $ = <I>(A;x) solves the 
following RH problem for x > 0. 

RH problem for <I> 



^»(2j-l) . 

shown in Figure [3j 



(a) <I> : C \ U|^;^e * — )• C is analytic, with the rays e 3 M+ oriented as 



7ri(2j-l) 

(b) $ has continuous boundary values on U|^]^e 3 \ {0}, and they are related 

by the jump conditions: 

$+(A) = $_(A) (^J G(A;x)-i^^ as A G M+ U e^M+, (3.11) 

$+(A) = $_(A)f , °V as Aee¥M+Ue¥R+. (3.12) 



-G(A;x) 1, 

(c) $ has the following behavior as A — t- oo: 

^{X) = I + 0{\-^). (3.13) 

(d) $ is bounded near 0. 
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Figure 3: The jump contour and jump matrices for <I>. 



We will prove the following results. 

Proposition 3.1 (i) //Re a > the RH problem for ^ is uniquely solvable for all 
but possibly a finite number of positive x-values {xi, . . . , Xk}, where xj = Xj{a, (3) 
and k = k{a, /3). 

(ii) If a > — ^ (Ima = 0) and Re/3 = 0, the RH problem for <1> is (uniquely) solvable 
for all X > 0. 

(iii) If He a > the asymptotic condition ^3.13^) for $ is valid uniformly for x G 
(0, +oo) provided that x remains bounded away from the set {xi, . . . , Xk}- 

Statements (i) and (ii) follow immediately from Theorem 11.81 (which will be proven in 
Section H]) . The third statement will follow from our asymptotic analysis of the RH 
problem for $ in Section HI 

We will now transform the jump matrices for <I> into the jump matrices for S near 1. 
Note first that the off-diagonal entries of the jump matrices for $ have branch points 
at ±1, and the ones for S at e^*. Let us therefore define a conformal mapping \{z) in 
a neighborhood of 1 which maps e~* to — |, e* to |, and 1 to 0: 

A(z) = Jln(z), zGf/. (3.14) 

Here we take the branch of the logarithm such that Inz > for z > 1, and the branch 
cut is along the negative real axis. We will furthermore need that e^^^-* = z"", and 
therefore set 

X = 2nt. (3.15) 



Let us choose the contours Ei and S2 near 1 in such a way that A maps Si U S2 
form 



onto the jump contour U^.e 4 M+ for We look for the parametrix P in the 



P{z) = E{z)<i>{X{z);2nt)Wiz), (3.16) 
where E is an analytic function in U, and W is given by 



1 , 



^„.(,) = -G(M.)n".."n^-.".„ for |.| < 1, 

\G{X{z)y2-^z2-^f{z)2-^ai, for|z|>l, ^ ^ 
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with cji ~ oj ■ Note that the branch points of G cancel the ones for / in U, and 

W is analytic in U \ C. 

If E is analytic in [/, it is easy to check using (I3.1ip - (j3.12j) that P{z) satisfies the 
same jump conditions as the matrix S with the jump matrices given in (j3.3p - p.5p . 
Since we evaluate ^{\;x) at x = 2nt, we need to impose the condition that 2nt does 
not belong to the set {xi, . . . ,Xk} of values at which the RH problem for <I> is not 
solvable. 

To fix E{z), let us consider the behavior of P on dU. From ()3.14p one observes that 
there exists c > such that for any < t < to 

|A(z)| > cn, z£ dU. (3.18) 

As n — )■ oo and if 2nt stays bounded away from the set {xi, . . . , Xk}, we can thus (by 
Proposition 13. ip use the asymptotic behavior (|3.13p for <I> to conclude that 

P{z) = E{z){l + 0{n-^))W{z), asn^oo, (3.19) 

uniformly for < t < to and z G dU. If to is sufficiently small, we can assume that 
e^* lie inside U and at a distance bounded from below away from dU. Then we obtain 
from (|3.17p and ()3.10p that (here and in ()3.22p below, 0(1) is a scalar matrix element) 



I Oil) \ 

z < 1 



W{z) = n 



Oil) 








Oil) 





Oil) 


0(1) 






(3.20) 

\z\ > 1 



as n — oo uniformly for < t < to and uniformly for z G dU \ C. 
Now set 

Eiz) = Niz)Wiz)~^. (3.21) 

One verifies directly, using the jumps for N and W across C, that E is analytic in a 
full neighborhood U oil. Furthermore, by (j3.20p . ([3^1 



as n — )• oo uniformly for < t < to and z € dU . 

Using this result and (|3.19p . we obtain the following matching condition on z G dU: 

Piz)Niz)-^ = Eiz) (/ + Oin-^)) Eiz)-^ =1 + n-^''^Oin-^)nl^"K (3.23) 

as n — )• oo uniformly for < t < to and z G dU . Note once again that the matching 
holds true if 2nt remains bounded away from the set {xi, . . . , Xk}- 

Since P has the same jumps as S inside U and 5(z)P(z)~^ = C'(ln(z — 1)) as z — 1 
for t > 0, it follows that the singularity is removable and Siz)Piz)~'^ is analytic in U . 
For later use, we note that 

Eie^) = n-'^''^g''^ai, £^(e"*) = n~^'^3/i'"3^i, (3.24) 
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Figure 4: The contour S/j. 



with g and h given by 

9 = (^) e-/^e- exp (-^ J] y.e*'^ + V,e'' ] , (3.25) 

:) ^ e-/^ exp (-1 V.e-^'^ + 1 y.e^^'^) . (3.26) 

^ V fc=-oo A:=0 / 



3.5 Final RH problem 

Define 

R(z) = { ^ ^ ^ \ ' ^ ' (3.27) 

[nl^''^S{z)P{z)-^n~l^''^, for z G [/. 

Using the RH properties of S, N, and P, we obtain the following. 

RH problem for R 

(a) R is analytic in C \ S/j, where T,r is the union of dU and the parts of Si, S2 
lying outside U (see Figure H]) . 

(b) R+{z) = R-{z)Jr{z) for z G Sr, where 

J^(z) = n^'^3p(z)^^(z)-^n-^'"^ for z G dU, (3.28) 

J^(^) = n^'^3A^(z)Jfc(z)iV(z)"^n~^'"3^ for z G Sfc outside C/, (3.29) 

and Jkiz), k = 1,2 are the jump matrices (j3.3p . (j3.5p of 5". 

(c) As z ^ 00, i?(z) = / + ©(z-i). 

Using (|3.3p . (j3.5p . and (j3.23p . we observe a crucial fact: 

Jj^{z) = 1 + 0{n-^), for z G (3.30) 

J^(2) = 1 + 0(6"'='^), for z G S/j \ aC/, (3.31) 
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as n — )■ oo uniformly in z and uniformly for < t < to ^-s long as 2nt remains bounded 
away from the set {xi, . . . , x^}. Thus the jump matrix J^^ tends to the identity matrix 
as n — )• oo. The RH problem for i? is therefore a so-called small-norm RH problem, 
and by a standard analysis (see, e.g., fl3i) we obtain the following statement. 

Proposition 3.2 Let < t < to. Then 

(i) The RH problem for R is solvable for n sufficiently large as long as 2nt remains 
bounded away from the set {xi, ... , x^}. 

(ii) If n ^ oo, 

R{z) = I + 0{n-^), uniformly for z eC\^R (3.32) 
and for < t < to such that 2nt remains bounded away from the set {xi, . . . , x^}. 

3.6 Asymptotics for ^lnD„ 

Reversing the transformations S ^ R, T ^ S, and Y ^ T hy (|3.ip . (j3.6p . and 
(I3.27p . we obtain the asymptotics for Y{z) under the conditions of Proposition 13.21 In 
particular, we obtain 

Y{z) = n-^''''{I + 0{n-^))n^''^P{z)z'"'\ for z near e*, (3.33) 

Y[z) = n-^^'-'il + C'(n"^))n^'"3P(z), for z near e"*. (3.34) 

as n — )• oo, uniformly for < t < to ii 2nt remains bounded away from {xi, . . . ,Xfc}. 
Using the definitions (13. 7p . (j3.16p . and (I3.24p . we find for P in the above formulas: 



P(z) = E{z)<^{z)W{z) = Dizf^ ( ° M W{z)-'^^{z)W{z), for z near e 



P(z) = E{z)^{z)W{z) = D{zy^W{z)-^^{z)W{z), for z near e\ (3.35) 

I 1^ 
1 0, 

(3.36) 

We will now substitute the asymptotics we obtained for Y into the differential 
identity (|2.8p for lnP)„(t). First, consider the case of z close to e*. By (|3.33p . we obtain 



Z 



(3.37) 



Using (j3.35p and (j3.17p . we further obtain 



P-^P', = -(73^ + W-^<^>-^<^>',W - ly-i^-Vg^VF + ^) , (3.38) 

where we defined AiyZ~) by the formula 

W{z) = A{z)''^ai. 

Expressions (j3.17p and (j3.8p give 

A' -t, a + B , a — B , f 1 e~* \ 1 ^ ^ , s 

^(e*) = ^e"* + —r^^~' 7 + -TT + (3.39) 



^ ' ' 4 4 \t sinh t / 2 

P) ^ 2sinht 



^ , y fcy,e('=-i)*. (3.40) 
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Therefore, we finally have for the 22 matrix element of P ^P^. at the point e*: 





^ ^ fc = -oo / 



(3.41) 

Now using the definition of it is easy to conclude that 



where is bounded in n as long as ) is bounded. Thus, we obtained the asymp- 
totic expression 

e\Y-^Y^h2{e') = -n + e\p-' P'^Uie') + (c|-^(t)0(l/n)l>(t))^^ (3.42) 

uniformly for < t < to as long as 2nt remains bounded away from the set {xi, . . . , x^}, 
with the second term on the r.h.s. given by (j3.4ip . 
Similar calculations at e~* give 

e-*(y-iy;)22(e-*) = e-*(p-ipi)22(e-*) + [^-\t)0{l/n)m) , (3.43) 

with 

e-\p-'P'Me-') = - ^ (7 + + \--'y^(^-') + e-(cl>-ici,;),,(e-) 



4 ^ 4 U 



^ fc=i / 

(3.44) 



Collecting (j3.42p and (|3.43|) together, substituting into (|2.8|) . and noting that 

^ z ^ 

gives 

Proposition 3.3 Let 

Hx) = -^{^-'<^'Mx/2) + ^($-i$l)22(-x/2). (3.45) 

T/ien 



dt = + 2^ 2^ U ^ sinht 7 

_^eV;(e*) + ^e-Vi(e-*) + 2nu;(x) 
Q + /3la + /3 a — e 



-0(l/n)$(x), 
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(3.46) 



where the error term is uniform for < t < tQ as long as 2nt remains bounded away 
from the set {xi, . . . , Xk} , and $(x) depends on ^{x/2), <!>(— x/2), a, j3 only and is 
hounded when these parameters are in a compact set. 

In the next section, we will analyze the Painleve functions <&(A). Namely, we will 
obtain their behavior at x = and oo which will be used in the last section to prove 
Theorems 11.11 and 11.41 



4 Model RH problem near 2=1 and the fifth Painleve 
equation 

Recall the RH problem for ^ stated in the introduction on the contour T given in 
Figure [J for j; > and Re a > — ^. In this section, we will analyze the ^-RH problem 
asymptotically for x near zero and infinity and prove Theorem 1 1 . 8 1 and Proposition [3Tj 
Moreover, we will give an explicit formula for w (defined by ()3.45p ) in terms of the 
Painleve V function v. 

For simplicity, we will omit the dependence of ^' on a, /?, and x in our notation when 
convenient. The behavior of ^ near and 1 can also be characterized in a different 
way from conditions (ll.SQh . (I1.40p . The following statement holds. 

Proposition 4.1 Let ^ = ^'(C; a, /3) satisfy the conditions (a), (h), (dO), and (dl) 
of the RH problem for ^ . Set 

*o(C) := ^(C)(C - l)^'^«r'^"^ (4.1) 

a-,8 

Then the function is analytic near and near 1. The branch cuts for C, 2 '^3 ^^^^ 
{C, — 1) 2 Q,j,g QfiQggfi Ji^re along [0, +00) and [l,+oo), respectively. 

Proof. The fact that \I'o is analytic near and 1 can be verified using the jump 
conditions for ^ for ((^ — 1) 2 '^^ and C, 2 . The isolated singularities at and 1 
are removable because of the conditions (dO) and (dl). □ 

Recall the function <I> defined in terms of ^ in (|3.9p - (|3.10|) . It satisfies the RH 
problem given in Section [3. 4[ 

We will now perform an asymptotic analysis of the RH problem for $ as x — )• +00 
and as x \ 0. 

4.1 Asymptotics for $ as x — t- +00 

Consider Figure [5] and define 

$(C;x) = $(xC;x), (4.2) 
for C, outside of the two triangular regions A and B. In these regions set 

$(C;x) = $(xC;x) G{xC,^) ^ ^ region B, 

$(C; x) = ^>(xC; x) {^Q(^^Q. ^-^ ' in region A. 

Now <I> is defined in such a way that it has its jumps only on the solid lines in Figure 
[5l We have 
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Figure 5: The jump contour for <^>. 

RH problem for $ 

(a) $ : C \ (fi U fa) ^ C^^^ ig analytic. 

(b) $ has continuous boundary values on Ti U T2 related by the conditions: 

$+(C) = $-(C)(j ^^"^^'^""^''y forCefi, (4.3) 

$.(C) = $-(C)(_^4^) ?)' ^'-'^ 

(c) $ has the following behavior as ^ — ^ 00, 

$(C) = / + 0(r^). (4.5) 

The jump matrices for $ are exponentially close to / as x — )• +00 because of the 
exponential factor in the definition of G (see ()3.10p ). Indeed, let us denote G'^(C) the 
jump matrix for the function ^{z), i.e. 

and 

Then, the following estimates hold: 

ll^-G|.llL.(f)' ll^-^$llL^(f)<^«"^^'"^^ 0<6<1, 00, (4.6) 

where 

r = fiur2, 

and we assume that the vertical parts of the contours Fi and F2 are given by the 
equations Re C = 1/2 — e/2 and Re^ = —1/2 + e/2, respectively. Estimates (|4.6|) imply 
that the RH problem for $ is a small-norm RH problem for large x and is therefore 
solvable in a standard way (see e.g., [13]) for x sufficiently large. Moreover, the solution 
<I> admits the integral representation. 
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with the function p{Q which is L2 - close to the identity, namely, 

l|/-p|li,(f)<Ce-f(i-). (4.8) 

The large x solvability of the <^>-problem means that the RH problems for <I> and ^ 
are solvable for x sufficiently large as well. In addition, we have that 

$(C;x) = / + ©(^e-f^^-^)), asx^+00. (4.9) 

This estimate holds uniformly for (" off the jump contour and this implies, in particular, 
that the following asymptotics hold as x — )• +00: 

$(±^;x) = H±^;x) = / + 0(e-§(i-^)), (4.10) 

cl>',(±|;x) = 0(e~f(i-)). (4.11) 



Furthermore, (14. 9p implies that 

$(A;x) =/ + C'(A"^), as A ^ 00, uniformly for X > C, C > 0. (4.12) 

The integral representation (j4.7p in conjunction with the estimate ()4.8p allows to 
evaluate the asymptotics for the Painleve function v{x) defined in (jl.4ip . Let Ci be 
the first coefficient in the large C expansion of the function 

$(C) = /+^ + 0(r'), C^oo. (4.13) 

(Note that this expansion is uniform for a; > C, C > 0.) Recalling the relation between 
the functions and ^, we obtain the following expression of the first coefficient Ci 

of the series (I1.38P in terms of the coefficient Ci , 

= ^i^ag + e-t'^«x^"«C'ix-^'^»et"^ (4.14) 
Together with (jl.4ip . this means that 

vix) = -Ci,u - xCi,i2Ci,2i. (4.15) 
On the other hand, from (j4.7p we obtain 

Ci = -J- Lp{C)(g^{C) - l) dC (4.16) 



2vri _ 

which leads to the estimate: 

1 , . \ 1 



^ [(G.{C)-l)d(: + 0{e-^('-^^), 0<6<^. (4.17) 



The estimate (j4.17p implies the following asymptotic representations for the entries of 
the matrix Ci as x — t- +00: 



Ci,i2 = -r^ / G(xC;x)-i(iC + 0(e-^(i-^)) 



-2/3. 



e 



-27ri/3 



-x--e-2 Ml -a-(3,2-2(3;x) + 0{e-^('-^^) 

T{a + p) 
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Figure 6: The jump contour and jump matrices for ^. 



T{a- p) 



?/^(l-a + /3,2 + 2/3;x) + C'(e- 



-x(l-e) 



-l+a+/3 



e 2 



r(a-/3) 



l + O ( - 



and 



Ci. 



11 



-Ci, 
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(4.19) 



(4.20) 



where ip{a,c;x) denotes the confluent hypergeometric function. The last estimate can 
be improved with the help of the differential identity ()4.96p which will be proven in 
Section I4.3i Indeed, this identity implies that 



ax 



Cl,12Cl,21, 



and hence 

Ci.ii 



^-2+2Q:g-a; 



T{a- /3)T{a + /3) 



1 + Oi- 

X 



(4.21) 



Substituting the estimates (|4.18|) . (|4.19|) . and (|4.2ip into the formula (|4.15|) . we arrive 
at the following asymptotic equation for the Painleve function v{x): 



v{x) = X 



-l+2a-x 



-1 



r(a - /3)r(a + /3) 



l + O 



(4.22) 



4.2 Asymptotics for and $ as x \ 



Write ^'i, . . . , for the analytic continuation of ^ from the indicated in Figure [T] 
sectors I, . . ., V to C \ [0, +oo), respectively, and consider the function 



^'(A;x) := e^'^^x 



X < 



^-11(^1;^)' 

^iii(| + l;x), 
^'iv(| + l;x), 
[*v(^ + l;x). 



for A in region I', 
for A in region II', 
for A in region III', 
for A in region IV', 
for A in region V, 



(4.23) 
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where the (modified) regions I', . . ., V are indicated in FigureEl so that ^ = ^{X; x, a, /3) 
has jumps on a contour which is partially shifted compared to the one for ^{X/x + 1). 
In particular, the intersection of the contour lines is now at A = instead of A = —x/2. 
From the RH problem for ^, one easily derives the RH conditions for ^. 

RH problem for 

(a) $ : C \ (^e^'-TR u M+ U [-x, 0]j C^''^ is analytic. 

(b) ^ has continuous boundary values on e^~M U U [— 2;,0] \ {—x,0}, and they 
are related as follows (with the orientation of the contour as in Figure [6]), 

as A G e-fM+, (4.24) 
as A G e^R+, (4.25) 
as A G eT^M+, (4.26) 

as A G eT^M+, (4.27) 

as A G M+, (4.28) 
asAG(-x,0). (4.29) 

(4.30) 



(4.31) 



(4.32) 

As A — )• in the other sectors, the behaviour of ^ is obtained by applying the 

^ \a + l3\ 

jump conditions to (14.321) . ^'(A) = 0{X 2 ) always holds. 

For small values of x, we will now construct a global parametrix and a local parametrix 
near for ^ and match them on the boundary of an e-neighborhood of A = 0. These 
constructions will lead to the uniform asymptotics for ^ and ^' as a; \ 0. 

4.2.1 Construction of the global parametrix 

Consider a fixed e-neighborhood of X = containing, in particular, the [—x, 0] part 
of the contour. Outside of this neighborhood, we expect to model ^ by the global 
parametrix M = M(A; a, j3) independent of x and solving the following RH problem. 



^>+{X) = ^^{X)[^^ ^ 
$+(A) = $_(A)(_^J(„_^) \ 
$+(A) = $_(A)(^^,(L/3) J), 

^+{X) = ^!^{X)[^^ ^ 

$+(A) = $_(A)e2^^^'"^ 
$+(A) = $_(A)e-^^("-^)'"3^ 

(c) As A — )• 00, 

$(A) = (/ + 0{X-^)) A-'^'^Sg-iAaa ^ 

(do) As A ^ -X, 

A + X M~ A + x 2" 



a-P a-13 
A + X 2 A + X 2 



\^(A) = O 
(dl) As A — )• in sectors V and V, 

\|A| 2 |A| 2 
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RH problem for M 

(a) M : C \ (^e=^xMu M+) ^ C^^^ analytic, 

(b) M has continuous boundary values on e^^MUM^ \{0} related by the conditions 

(1 „7ri(«— 

^ ^ J , as A G eTM+, (4.33) 

M+(A) = M_(A)(_^_i(,_^) J), asAGe^M+, (4.34) 

M+(A) = M_(A)(^^,(L;3) i), asAGe^]R+, (4.35) 

Q j , as A G e^R+, (4.36) 

M+{\) = M^{\)e^^''^''-\ asAGM+. (4.37) 

(c) We have 

M(A) = (I + C'(A-^)) A-^'^^e-^^^s^ as A ^ oo. (4.38) 

We can solve this RH problem explicitly in terms of the confluent hyper geometric 
function. Inspired by the constructions of [23l [H] , we define 



g-i7r(2/3+a) \ _hl / A°V(a + /3, 1 + 2a, A)e*'^(2/3+a) 

L ^ 



^^^^ V e-(/^+2") j ' l^A-V'(l - a + /3, 1 - 2a, A)e-(^-3") r^^+^f^ 

A-V(l + a - /3, 1 + 2a, e-A)e-(/^+-)qi±f^^ .^.Sg-A^s/s 
A-^V'l-a - /5, 1 - 2a, e-*'^A)e-*™ 



(4.39) 



where 'il^{a, b, x) is the confluent hypergeometric function, and T{x) is Euler's F-function. 
Furthermore let 

M(A) = Mi(A):=i7(A)f J, forO<argA<f, 

M(A) = M2(A) := H{X), for f < arg A < ^, 

M(A) = M3(A) := H{X) (^_^.X-l3) i) , for f < arg A < f , 

M(A) = M4(A) := //(A) ^^^^^J^^ _ , for ^ < arg A < ^, 

M(A) = M5(A) := F(A) 

1 _e-'r*("-/3) \ .^^ 

2isin7r(a-/3) -2ie--(-« sin7r(a - /?) + 1 j ' f < ^rg A < 27r. 

Using known properties of the confluent hypergeometric function, one verifies as in \lT\ 
that M satisfies the prescribed RH conditions. 
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In order to match M later on with the local parametrix near zero that we construct 
in the next section, we will now need to rewrite M in the form in which the structure 
of its singularity at A = becomes more apparent. 

Recall the following properties of the confluent hyper geometric function (see, e.g., 

6(a r z) p*™ „-iiT{c-a) 

' ' ^ - i'{a,c,z) + \, , tP{c-a,c,e~'^z)e', (4.40) 



r(c) r(c-a)^' ' ' ' T{a) 
i^ia, c, z) = z'^'^i^ia - c + 1, 2 - c, z), (4.41) 

i^ia, c, z) = — rTT'/='(a) c, z) + ^K ^'^'ifia - c + 1,2 - c,z), c ^ Z, 
i (a — c + 1) L [a) 

(4.42) 

^p{a,c,z)=e^ip{c — a,c,—z). (4.43) 

where 

/ ^ 1 a(a + 1) • • • (a + n - 1) / ri -, n ^ a a a\ 

n=l 

is an entire function. Let us focus on the region III'. Assume first that 2a 7^ 0, 1, 2, . . . 
(noninteger 2a combined with our general condition Re a > —1/2). Using the proper- 
ties (|4.40p and (|4.4ip for simplification of the first column of M3, and (|4.42p . (|4.43p for 
the second, we easily obtain 

M3(A) = ^(A)A°'^^G3, 2a/ 0,1,..., a ± /3 / -1, -2, . . . , (4.45) 

with the branch of A^° chosen with < arg A < 27r. Here 

p,,, _ ,-A/2 ( --''^^^'^ + /3, 1 + 2a, A) 

" \-e-^-'^m^^i^ + a + 1 + 2a, A) 



r(l+2a) 

e*"^""''^rT^</'(-« + Al-2a,A) ^ 
e^"("+^)T^<^(l-« + /3,l-2a,A)^ 



(4.46) 



is entire, and G3 is the constant matrix 



Applying to (j4.45p the jump conditions we readily obtain the general formulas 

M,{\) = E{\)X'"^Gj, j = l,...,5, 2a/0,l,..., a ± /3 / -1, -2, . . . , 

(4.48) 

with appropriate constant matrices Gj. In particular, G3 is given by (j4.47p and 

A + cie-^^("-^) _eW"-/3)\ 
^i=( e-W-/3) )■ (^-^^^ 

Consider now the case when 2a is an integer. In this case, the calculations for the 
first column of M3 remain the same, whereas for the second column (note that (j4.42p 
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does not hold now) we use the known logarithmic formulas for the -(/'-function. We then 
obtain 



M,(A) = i?(A)A"'^3 (^J i = l,...,5, (4.50) 

7(A) = ^ sin7r(a + /3)ln(Ae~*^), if 2a = 0,1,..., (4.51) 

vr 

where E(X) is analytic at zero, and Gj are constant matrices. In particular, G3 = I. 
4.2.2 Construction of the local parametrix near 

We construct a local parametrix P in Ug in such a way that it has the singularities and 
jumps of ^ and matches with M to the main order in x on the boundary. We have 
to consider the cases 2a ^ Z and 2a € Z separately because of the different behavior 
of M at dUir. As examination of the final formulas show (see (j4.6ip . (j4.64p . ()4.7ip . 
and (j4.72p below), our constructions in this section will be valid for all a, /3 such that 
a lb /? 7^ —1, —2, . . . This is exactly the restriction on M in the previous section. Some 
preliminary expressions, however, are valid under stronger conditions a ± /3 ^ Z (cf. 
()4.56p ). We do not mention these conditions as they disappear in the final formulas, 
namely: the singularities in a, (3 of J and J defined below cancel with the zeros of cq, 

We first deal with the case when 2a 7^ 0, 1, 2, ... . Since we have a problem with 
2 singular points with power-law behaviour at —x and and require a power-law be- 
haviour at the boundary dUg (see (I4.48P ) we expect from the general principles (e.g., 
|16j ) a parametrix in terms of the hypergeometric function. Indeed such a parametrix 
was found by Jimbo [25] in the generic case of Painleve V equation. Instead of trying 
to specialize it to our situation, we provide a direct construction below. For 2a ^ Z , 
define P = P{X;x,a, 13) in 17^ by the expressions 

P(A) = Pj{X), with j = 1 for A in sector I', j = 2 in sector 11', and so on, (4.52) 

where 

P,{X)=E{X)(l "°'^(^=;^'«'^))(A + x)^'^3;,^.3G^. (4.53) 

with the argument of the roots between and 2tt. Here E is given by ()4.46p . Gj are as 
in Km . 

J(A; X, a, (3) = - + " 1^1^ " ~ l+a+f3, 2+2a, e^x/A), (4.54) 

vr —A i (2 + 2a) 

CO = -e^-^"^"^" + ^^^^""^""^\ (4.55) 
sin 27ra 

and F{a, b, c, z) is the hypergeometric function of z with parameters a, b, c. The 

argument of z is chosen between and 2tt. For c 7^ 0,-1,-2,..., this function is 
represented by the standard series 

, , _ ^ , ^ a(a + 1) • • • (g + ra - 1)6(6 + 1) • • • (6 + n - 1) 

r [a, 0, c, zj — l+> r , i\ I' 

^-^ c(c + 1) • • • (c + n — 1) n! 

71 = 1 
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converging in the disk \z\ < r < 1 of any radius r < 1, and is extended to the analytic 
function in the plane with a cut [l,+oo). Therefore, the function F{a,h,c,e™ x/ \), in 
particular the one in (j4.54p . is analytic in A-plane outside [— x, 0]. We will now find the 
jump of J(A) on [— x, 0] and the structure of its singularities at —x and 0. First, using 
the transformation of the hypergeometric functions from the one with the argument z 
to those with the argument 1/z, we can write: 

F(l, 1 + a + /3, 2 + 2a, z) = - . 1" ^ ^ ^1\V('^\ m 

sm7r(a + pj r(l + a + p)T(l + a — p) 

X (e'^z-^) F(l + a + /3, -a + /3, 1 + Q + /3, 1/z) (4.56) 

+^-^e'^z-^F{l, -2a, 1 - a - /3, 1/z). 
a + p 



Since the function 

F(l + a + /3, -a + 1 + a + /?, 1/z) = (^1 

and F(l, —2a, 1 — a — /3,1/z) are analytic outside [0, 1], we have on (1, +oo): 
F(l, 1 + a + /3, 2 + 2a, z)+ - F(l, 1 + a + /3, 2 + 2a, z)_ 

_|^|-l-Q-/3 



1 

1 - - 

Z 



r(l + a + /3)r(l + a - /3) ' ' 
Setting here z = e^'^x/X and using (|4.54p . we obtain 

J+(A) = J_(A) + 2i|Ar+'^|A + x|"-^, Ag(-x,0). (4.57) 
By ()4.48p we can write 
P(A) = ^(A) (^J ""oJiX)^ ^ ^^f^,3^_ii_£,3^^^^_i^^^^^_ ^^_^g^ 

Since the product of the factors to the left of M is analytic in C\ [—x, 0], the expression 
()4.58p implies directly that P satisfies the jump conditions for M, see (|4.33p ~ (j4.37p . 
except on (—x,0) where M is analytic but P is not. For A G (— x,0), it is convenient 
to use (|i33]) . (07D . and (|i371) to verify that 

P+(A) = P_(A)e-^^("-^)'"3. 

This is only true if Cq is given by (j4.55p . We have thus constructed P in such a way 
that it has exactly the jump conditions for ^. 

From the fact that P and ^' have the same jumps, it follows that "^P^^ is analytic in 
Ue except possibly at the points —x and 0. Let us investigate the behavior of near 

a-f) a + fi 

these points in some detail. Recall that by Proposition |1]T] ^'(QC 2 '^3(^i^_i^ 2 is 
analytic at and 1, which implies by (|i:23l) that $o(A) := $(A)(A + x^^'^'^X^''^ 
is analytic at —x and 0, where at —x the variable A is in the region ///', while at 0, 
the variable A is in the region /' or V' (for other regions at 0, the appropriate jump 
conditions should be applied to ^(A)). Consider first A close to 0. By (14.53p . we have 

$(A)P(A)-^ = $o(A)(A + x)'^'''X-^''''G^\X + x)-^'"'-' X'^'"'' 



1 —coJ{X;x,a, (3) 




i""'"'^^^ E(A)-\ AG/'. (4.59) 
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Substituting (j4.56p with z = e^'^x/X into (|4.54|) . we obtain the following representation 

J{X;x,a,(3) = . \ x^"^ ^ ' F{1, -2a,l-a- (3, -X x). 

sin7r(a + pj vri (1 + 2a j 

(4.60) 

Therefore, 

CO J(A; X, a, (3) = -e^"" "^'^^^" ~ A°+^(A + x)""^ + Fo(A), 

sin 27ra 

7? ^2a,2^ia Sin7r(a - ^) r(l + g - /3)r(-2a) « \ / ^ 

Fo(A) = -X e — F(l,-2a,l -a-/3,-A/x). 

TT i (i — a — pj 

(4.61) 

Since -Fo(A) is analytic close to A = 0, this representation gives an explicit expression for 
the singularity of J(A) at 0. (Note that the arguments of A, A + x here are between — tt 
and vr as they originate from the hypergeometric function above, whereas the arguments 
of A, A + X in the first row of (j4.59p are between and 27r.) Substituting ()4.6ip into 
(|4.59p we see that the singularity cancels, and ^'(A)P(A)~^ is analytic at A = 0. 

In order to analyze the singularity of J at A = — x, apply first the transformation of 
the hypergeometric function between arguments z and 1 — z to ()4.56p with z = xjX. 
We then obtain: 

F(l, 1 + a + /?, 2 + 2a, -x/A) = 

7re-^"" r(2 + 2a) / A\ ^+"+^ / A\ ""^ 

sin7r(a — /3) r(l + a + /3)r(l + a — /3)\x/ \^x 

) (4.62) 

-^^^-F(l, -2a, 1 - a + /3, 1 + A/x), 
a — p X 

and therefore by (j4.54p . 

e-2™A°+''(A + x)°-^ 



J(A; X, a, /3) 



sin7r(a — /3) 



(4.63) 



+x^" ^^";g|f;+^; + ^^ F(l, -2a, 1 - a + /3, 1 + A/x). 
Using ()4.55p . we obtain 

CO J(A; X, a, /?) = ^^'^^^^ + ^^ °+^(A + x)""^ + Fi(A), 
sm 27ra 

T?(\\ 2a 2^ia Sin7r(a + /3) r(l + a + ^)r(-2a) 

Fi(A) = x e — — r F(l,-2a, 1 -a + /3, 1 + A/x). 

TT i (1 — a + pj 

(4.64) 

This representation explicitly displays the singularity at — x because F\{X) is analytic 
near A = — x. In the same way as for A = 0, we now obtain that ^'(A)P(A)~^ is analytic 
at A = —X as well. 

For A at a fixed distance away of the origin, say for A G 9?7e, it follows from ()4.58p 
that 

P(A;x)M(A)-i =I + C)(x) + 0(xi+2a)^ as x \ 0. (4.65) 
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We now consider the case 2q! G Z. We again set 

P(A) = -Pj(A), with j = 1 for A in sector I', j = 2 in sector II', and so on, (4.66) 
but now with 

P.(A) = ^(A)(; ^^f))(A + x)^-A^-(; (4.67) 

where 



C2 = sin7r(a + /3) sin7r(Q; — /3), 

TT 



(4.68) 



J(A; X, a, 13) 



with J(A) given by (I4.54p . 
Note that 

J^{\) = J_(A) + 2i|Ar+^|A + a;|"-'^ln|A|, A G (-x,0). 



(4.69) 



(4.70) 



A similar derivation to the one above shows that P again satisfies the same jump 
conditions as the ones for To analyze the structure of the singularities of P(A) at 
and at —x, we need, as above, to find suitable expressions for C2J(A) at these points. 
Applying the differential operator ^{d/da + d/dfi) to (|4.60p (it is convenient to write 
)a+p ^ gi7r(a+/3)(g-i7r^-)a+/3 ^^.g^^ ^YiQn taking 2a G Z and using the fact that 



sin7r(a — /3) 
sin7r(Q + /3) 

we obtain 



(-1) 



2a+l 



2a G Z, /3 ^ 



C2J(A;x,a,/3) = 
C2 ( d 



(-1)2°+! (Ae-'") 



i-!T\a+f} 



TT 



sin^(a-/3)A"+'' ln(Ae- 



(A + x)"-^ + Fo(A) 



X 



2a 



T{a + f3)T{l + a- (3) 
7rr(l + 2a) 



F(l,-2a,l -a-/3,-X/x) 



(4.71) 



As above, this expression can be used to show that ^'(A)P(A)~^ has no singularity at 
A = 0. To analyze a neighborhood of A = —x, note that the first term on the r.h.s. of 
(j4.63p can be written as 



-7ri(a— /3) 



sin7r(a — /3) 



(e-^"A)°+^(A + x) 



a-/3 



and application of the operator ^{d/da + d/d(3) to the fraction in the above formula 
gives zero by antisymmetry in a and /3. We finally obtain from (I4.63P 



C2J(A;x,a,/3) = - sin7r(a + /3)e-2^"A°+''(A + x)"-'' ln(Ae-*'') + Fi(A) 



\a-l3 - 



vr 



„2a 



r(a - /3)r(l + a + /?) 
7rr(l + 2a) 



F(l,-2a, 1 - a + (3,1 + X/x) 
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(4.72) 



which can be used to see that ^(A)P(A) ^ has no singularity at A = —x. 
It is easy to see that 



P{X;x)M{X)-^ = I + 0{xlnx), foiXedU, (4.73) 
as X \ 0. 

Note that using an integral representation for the hyper geometric function, we can, 
if Re (a it /?) > —1, represent J(A) and J(A) in the following form, which makes the 
jump conditions (|4.57p and (|4.70p obvious: 

T^J-x C-A vry_^ ^-A (4.74) 

Re(a±/3) > -1. 

Now define 

|$(A)M(A)-, forAeCXf/., 

^ ^ \^'(A)P(A)-i, forAGC/,. ^ ^ 

This function satisfies the following problem. 

RH problem for R 

(a) R is analytic in C \ dUe. 

(b) The jump condition for R is 

i?+(A) = i?_(A)(/ + e(A)), for A G dU, with e(A) = o(l) as x \ 0. (4.76) 

(c) R{X) = 1 + 0{X-^) as A ^ oo. 

For X sufficiently small, this is a small-norm RH problem, and it follows that the RH 
problem for R is solvable, say for < x < 6. Prom the invertible transformations 
I— I— ii and ^ i— <I>, it follows that the RH problems for and ^> are solvable as 
well for < X < (5. We also have R{X) = I + o(l) uniformly for A € C \ dU as x \ 0. 
In particular this holds at infinity, which means that 

i?(A) = / + O , asA^oo, (4.77) 

uniformly for small x. Tracing back the transformations ^' i— t- ^ i— ?• i? and i— )• we 
can conclude that 

$(A; x) = I + 0{X~^), as A oo, uniformly for < x < 5. (4.78) 
Moreover, using ()4.6ip we obtain for 2a ^ Z 

$(|;x)a3$(|;x)-i = e--'^^/^R{0)E{0) ("^^ '^^f^^^ E{Q)-^ R{Q)-^e-^-^'\ (4.79) 
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where Fq{X) is defined in (|4.61|) . Similarly, 

(4.80) 



where -Fi(A) is defined in (I4.64|) . 
Therefore, we have 

^.(±|;x)a3<i>(±|;x)"i = 0(l) + 0(x2"), x\0, 2a ^ Z. (4.81) 

Similarly, by (lOTD . KTTh . K72\i . 

$(±|;x)a3«>(±|;x)-^ = 0(l) + 0(a;2") + 0(x=^"lna;), x \ 0, 2a G Z. (4.82) 

We will now estimate w{x) given by p.45p as x \ 0. First, using the connection 
(|3.9p between $ and ^, we obtain 

. _ H)! (^ _ ^ , ,) . , . (4.3) 

Expressing = and using (|4.75p . and the formulas for P{\), we obtain 

after straightforward calculations that 



($-^$'Jn(x/2) = ^ + (di + (i2X^")(l + 0(x)), di 



X 



2a ' 



_ g - /3 r(l + g + /3)r(l + a - /3) r(-2a) 
^ 1 + 2g r(l - g + /3)r(l - g - /3) r(l + 2g)2 ' ' 



(4.84) 



(4.85) 



Similarly, 

(^-'^a)22(-x/2) = + (di + d2x2")(l + 0(x)), di = 

a + /3 r(l + a + /3)r(l + a-/3) r(-2a) 
^ l + 2ar(l-a + /3)r(l-a-/3)r(l + 2a)2' ' 

Substituting these expressions into (j3.45p . we obtain as x \ 0: 

, , a2 - /32 a2 _ ^2 

r ,, r(l + a + /3)r(l + a-/3) r(l-2a) 1 ] 

r(l-a + /3)r(l-a-/3)r(l + 2a)2l + 2a/^^ + ^^'^^^' ^ 

(4.86) 

Similarly, we verify using (g^ZD, ([iTTT]) . (ji:72]) that as x \ 0, 



w{x) = — — ^ + 0(l) + 0(x2°) + 0(x2"lnx), 2a gZ. (4.87) 
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4.3 Differential system for \1/ 

So far we know that there exist 6,M > such that the RH problems for ^ and $ are 
solvable for x > M and for < x < 6. We will derive differential equations for ^ with 
respect to x and C- This will lead to the Painleve V equation and will help us to find 
an identity for the function w given by ()3.45p in terms of v. Here we follow similar lines 
as in |16j Section 5.4]. 

From the RH conditions for ^, it follows that, for any x for which the RH problem 
is solvable, the (matrix) function A((!^;x) = /^{C', x)'$~^ {C; x) is a rational function in C, 
with simple poles at and 1. Indeed, A is meromorphic because ^ has constant jump 
matrices, A is bounded at infinity because of (jl.SSp . and has simple poles at and 1 
because of (j4.ip . Similarly, B{(;x) = "^^iC'i x)'^~'^{C; x) is a polynomial of degree 1 in 

It follows that ^ satisfies a linear differential system of the form 



^c(C;^) 



*(C;x), (4.88) 



^^iC x) = [i?i(x)C + Bo{x)] ^(C; x). (4.89) 

Substituting the large (^-expansion (jl.38p for into (j4.88p and (j4.89p . we can express 
the coefficient matrices ^oo, ^i; Bq, and Bi explicitly in terms of the entries of Ci 
and C2- 

Aoo = -^T3, (4.90) 

, _( —f3 + q + xrt —2f3r — xh + xrq + r + xr\ (a qt^ 

° " \ 2l3t + xj -xtq + t-xt (i-q-xrt J ' ^ 

_ f —Q — xrt 2f3r + xh — xrq — r\ ,^ , 

^ ~ \-2/3t - xj + xtq - 1 q + xrt ) ' ^ ■ ) 

Bi = -1^73, (4.93) 
Bo = (_°, , (4.94) 
where q = q{x), r = r{x), t = t{x), h = h{x), and j = j{x) are given by 

- (t1 -12)) • - 0(1-) <^-^=' 

(note that the trace of Ci must be zero since the determinant of ^ is equal to 1). 
Equating the C'(l/C)-terms in "i^^ = (BiC + -Bq)^ gives the identities 

q'{x) = r{x)t{x), (4.96) 
h(x) = -r'{x) + r{x)q{x), (4.97) 
j[x) = t'{x) + t{x)q{x). (4.98) 

Furthermore, equating the mixed derivatives = leads to the compatibility 
condition 

Ax - Bc + [A, B] = 0, [A, B]= AB- BA. (4.99) 
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Let us follow [TB] and write 

^(^) — — ~ ^^^^ ~ ^I'i^) ~ — ~ ~ xr{x)t{x), (4.100) 

xt 

u(x) = l + — -— — . (4.102) 

^ ^ (2/3 + 1 - x)t(x) + xt'(x) ^ ^ 

Using Proposition 14. 1 1 one shows as in |IB] that det Aq = — if^i^L and det Ai = — '^"^^^ . 
It then follows that the matrices Aq, Ai, and Bq can be written in the form (the elements 
(11), (22), and (21) of Aq and Ai are easy to verify directly, and (12) follows from the 
expression for the determinant): 



A„=h+!^ "^<"-°_J«), (4.103) 
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Writing the compatibility condition (j4.99p in terms of the functions u,v,y, one verifies 
that u, V, and y solve the system of ODEs 

xu^ = xu- 2v{u - if + {u- l)[(a - p)u - (5 - a], (4.106) 

xvx = uv[v - a + l3]--{v - (5 - a), (4.107) 

u 

xyx = y{^-2v + a + l3 + u[v-a + (3] + ^-x'^, (4.108) 

which is part of the content of Theorem 11.81 (iii). Eliminating v from the first two 
equations, one shows that u solves the Painleve V equation (jl.22p - (jl.23p . 
Define 

a{x) = xq{x) - '^-^x. (4.109) 

It follows from ()4.100p that 

a' = -V, (4.110) 
and therefore, by (|4.107p . 

-xa" = uv{v -a + /3)--{v- /3-a). (4.111) 

u 

Moreover, in view of (j4.96p . we have that 

a - xa' = -x'^q' = -x'^rt = a;^(i?o)i2(-Bo)2i- 
This equation can be rewritten with the help of (I4.105P as 

a — xa' = —^v — a — P — u{v — a + /3)^ ^ 

= uv{v - a + p) + -(v - 15 - a) -2v'^ + 2av. (4.112) 
u 
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Using (|4.11Up . (|4.11ip . and (j4.112|) . we can check directly that the function cr(x) satisfies 
the (T-form of the fifth Painleve equation ()1.2ip . 

The system (j4.88p - (|4.89p is the Lax pair associated with Painleve V. Since the RH 
problem for ^(<^;x,a,/3) is solvable for < x < 6 and for x > M, the Lax matrices 
AQ(x;a,P), Ai{x;a, f3), and BQ{x;a, /3) exist for those values of x. However, the system 
()4.106p - ()4.108p has solutions which are meromorphic in C \ {0} with a cut from zero to 
infinity, which implies that Aq, Ai, and Bq exist for all but (possibly) a finite number 
of positive x-values. Using appropriately normalized solutions to (|4.88p - (|4.89p . the RH 
solution ^ can also be constructed for all but possibly a finite number of positive x- 
values [16j. This proves Theorem ll.8l (i) and the equivalent statement for <I>, Proposition 
13.11 (i). Furthermore, the differentiability of ^ with respect to x, see (j4.89p . implies 
that the asymptotic condition ()1.38p . and thus also ()4.38p . holds uniformly as long as 
5 < X < M ii X remains bounded away from the set of x- values for which the RH 
problem is not solvable. Together with (j4.12p and (j4.78p . this proves Proposition 13.11 
(iii). 

Remark 4.2 The functions u, v, y appearing in ()4.103p - ()4.105p are particular solutions 
to the system (j4.106p - (|4.108p . Other solutions can be obtained by considering RH 
problems for ^ with modified jump matrices and modified behavior near and 1, 
corresponding to different monodromy data, see [2|. 

Remark 4.3 The RH problem for ^ is not the standard RH problem related to the 
fifth Painleve equation. In |16l I18j . a RH problem was posed on a contour UqUUiL) T, 
where Uq and Ui are small circles surrounding 0, and 1, and where r = M\ ([/oU?7i). 
The equivalence of a particular case of this RH problem with ours can be verified 
directly using Proposition 14.11 In order to avoid confusion with the notations in [16J, 
we note that the system ()4.106p - (|4.108p is written with parameters ^Oi ^ij and 6*00 in 
|16j . which in our setting are given by 



eo = -f3-a, 01 = a- /3, 0oo = 2/3. (4.113) 
Proposition 4.4 Set 

a{C;x) = {^{C;x)a3^-\C;x))^^. (4.114) 
Then the identities 

-^a{0; x) = Ao,ii = -v{x) + (4.115) 

^a(l; x) = -Ai,u = -v{x) + (4.116) 
hold, with V defined as before by ^Ljl^- 

Proof. Substituting ^ expressed from (j4.ip into the differential equation 

and comparing the residue of the left- and right-hand side at leads to an expression 
for ^'(C)<T3^'(C)~^ as C 0, in terms of Aq. By (l4.1U3p . this gives the first identity. 
Comparing the residues at 1 gives the second identity. □ 
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Proposition 4.5 Let w he defined by ^il.45\ )- Then 

v{x) = -{xw{x)y, (4.117) 
a{x) = xw{x), (4.118) 

r+oo 

a{x) = / v{C)dC. (4.119) 



Proof. It follows from Proposition 14. II that ^ can be written in the form 

^iC) = EiOC'^'^', ^(C) = F{C)iC - 1)~'^'^«, (4.120) 

with E analytic near and F analytic near 1. Let us write 

EiC) = EoiI + EiC + OiC^)), asC^O, (4.121) 

F(C) = Fo(/ + Fi(C - 1) + one - If)), as C ^ 1. (4.122) 

Substituting (j4.120p and (I4.12ip - (l4.122p into Kmh . we obtain the identities 

K,x = BoEo, E[,^ = Eq^BiEo, 

H,x = EqFo, F{,^ = Fq'^BiFo, 

which imply by Proposition 14.41 in particular, that 

E'lM^) = {E^^BiE^)22 = {E^BiE^%2 = ^«(0;x), Fl^^x) = -^a(l;x). 

(4.123) 

On the other hand recalling equation (j4.83p . we obtain 

*-'(-f-«(-f^^'). = -^^ + ^^- <^-'^^' 

*-'(f;.«(|;.))„ = i-^ + iF,„. (4.125) 

From (f05]) . it follows that 

and by (j4.123p together with Proposition 14.41 we obtain —{xw{x))' = v{x). From 
()4.110p and (|4.117p . it follows that (j{x) = xw{x) + constant, where a is defined by 
(|4.109p . To determine the constant, note first that, as follows from (j4.14p and (|4.17p . 
q{x) = Ci^ii — )• (a + /3)/2 as X — )• +00, and hence, a{x) — )• as x — )• +00. On the 
other hand, as follows from ()4.10p . ()4.1ip . we have xw{x) — t- as x — >■ +00. Hence 
the constant in question is zero, and we obtain (14. USD . Equation (I4.119P is obtained 
similarly. □ 

Combining (j4.117p . (jilM]) . (gSZl), and (j02]) . we obtain pliG]) . The expressions 
(OHD . gSZD, diHH]) . and (lilT9D imply IfLiT]) . 

Proposition 4.6 Let v he defined by ^1.41\ )- Then v(x) is real for x > i/Ima = 0, 
a > and Re/3 = 0. 
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Proof. Suppose that a > — ^, Re/3 = 0, and that is a solution to the RH 

problem for ^ given in Section [L3l Then it is straightforward to verify that the function 
^ defined by 

$(C) := ai^(-(C- l/2y)aie±^^^'^3^ if ±ImC > 0, 

with fji ~ ' solves the RH problem for '^{( + 1/2) for real x up to a constant 

factor. Therefore, by uniqueness, 

^'(C + 1/2) = C(x)'^^$(C), (4.127) 

where C{x) is independent of (. By (I4.114p . it follows that a{0;x) = a{l;x). Sub- 
tracting the complex conjugate of (|4.116p from (|4.115p . we conclude that v{x) = v{x). 

a 

4.4 Solvability of the RH problem for \E' 

In this section, we will prove Theorem 11.81 (ii) and Proposition 13.11 (ii): we will prove 
that the RH problem for ^ is solvable for all positive values of x if Re /3 = and 
Ima = 0, a > — ^. 



4.4.1 Vanishing lemma for Painleve V 

For a general class of RH problems, it is known that solvability of a RH problem is 
equivalent to the triviality of a homogeneous version of the RH problem |16t [2H] . For 
the case of Painleve V, this has been used in |T8] for a slightly different but equivalent 
RH problem (cf. Remark l4.3p . In our case a sufficient (and necessary) condition to prove 
the solvability of the RH problem for ^ is given by the following so-called vanishing 
lemma. 

Lemma 4.7 (Vanishing lemma for Painleve V) Let x > 0, Ima = 0, Re/3 = 0, 
and suppose that satisfies the RH conditions (a), (h), (dO), and (dl) of the RH 
problem for ^ , with condition (c) replaced by the homogeneous asymptotic condition 

^Q{C)e^^''^ = 0{C^), as C^oo. (4.128) 

Then it follows that = 0. 

Remark 4.8 A vanishing lemma was proven in [18] for a family of solutions to the 
system ()4.106p - ()4.108p . Our solution, however, is not contained in this family, and the 
vanishing lemma requires a different proof in our case. For the proof of the vanishing 
lemma, we follow similar lines as in [12, Section 5.3]. 



Proof of Lemma 14.71 Suppose we have a solution to the homogeneous RH 
problem. We will then prove that = 0. Let us first define a function M as follows, 

M(C) = ^o(C + ^)e^«"^ iff<argC<f, 
MiO = ^oiC + i)e-"^^"3efC-3^ for < argC < f , 

M(C) = ^'o(C + ^)e''*^'^3ei^'"^ for ^ < argC < 2tt. 
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for 


f <argC< f, 




for 


f <argC<^, 




for 


f <argC< f 


j g7ri/3(T3gf Co-3^ 


for 


f <argC< ^ 



In the remaining regions where j < argC < ^ or ^ < aig( < ^ we define M as 
the analytic continuation of M from the other sectors in such a way that M has jumps 
only on the imaginary and on a part of the real axis: 

1 /l _e'r«(a-^)\ 

M(C) = ^o(C+2)(; 1 ) 

M(C) = *0(C + ^)(_^J(a-/3) 1 
1 /l _p-'ri(a-/3)\ 

MK) = *.« + ^)(„ ^ ) 

Then M satisfies the following RH conditions. 

RH problem for M 

(a) M is analytic in C \ (iM U [-i, i]). 

(b) M satisfies the following jump conditions on the contour (iMU (—5, 5)), with iM. 
oriented upwards and (—5,5) oriented from left to right, 

M+(C) = M_(C)Fi(C), as C e (0,+ioo), (4.129) 

M+(C) = M_(C)F2(C), as C e (-ioo, 0), (4.130) 

M+(C) = M_(C)e--*("-^)'^3^ as ^ g (_i^0), (4.131) 

M+(C) = M_(C)e-'^^("+^)'^3^ as C e (0, i), (4.132) 

with 

J, (4.133) 

(c) For fixed x > 0, 

M(C) = ©((-'), as C ^ 00. (4.135) 



(dO) As C ^ -i 



2' 



^(0 = f(-||, (4.136, 



(dl) As C ^ +i 



2' 



ic-ii-t |C-i| 



M(0 = O(|J_||-| (4-137) 
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Let us now define a function -ff (C) in terms of M and its Hermitian conjugate as 
follows: 



H{C) = M(C)M*(-C). (4.138) 

Because of the condition (c) of the RH problem for M, we have that H{C,) = 0{C,^'^) 
as C — ^ CO- Furthermore, using the jump condition (j4.13ip - ()4.132p for M, we obtain 
that H has no jump across (0, ^). (This is only true if Re/3 = 0, Ima = 0.) Therefore, 
H is meromorphic for ReC > 0, with an isolated singularity at |, which is removable 
because of (|4.136|) and (j4.137p . Using Cauchy's theorem, we then have 

H^iC)dC = 0, / HliC)dC = 0. (4.139) 

o J —ioo 

Because of the jump conditions for M, the first integral implies that 

M^{C)Vi{C)M*_{C)dC + / M_{C)V{{C)M*_iC)dC = 0. (4.140) 

o JO 

Summing up this expression and the one obtained from the second integral in (j4.139p . 
we find, using (j4.133p . (|4.134p and the fact that x is real, 

|°^M_(C)(o 2e-;3)M:(C)dC + ^'"M_(C)(o ^e^) ^-(^^C = 0. (4.141) 

Since Re /3 = 0, it follows immediately that the second column of M_ is identically 
zero on iR \ {0}. From the jump conditions (I4.129p - (l4.130p . it then follows that the 
first column of M_|_ is zero on iR\ {0} as well. Therefore, we have that Mj2(C) = for 
ReC > 0, and Mji{() = for Re^ < 0. Let us now define 

yMji[Q, as ReC > 0, 

so that gj is analytic in C\ (fMU [— |, ^]). Furthermore, gj is bounded except near zt^. 
On iR, g has the following jump relation. 



93 AO = 93 AO X <^ _ ' ' '. ' (4.143) 

e e as C G (— ^oo,0). 



Now we write 'g for the analytic continuation of g from the left half plane to C\ [— ^ , +oo) , 

'g(C), asReC<0, 

?(C) = < c/(C)e™e-^^, as ReC > 0, ImC > 0, (4.144) 

,9(C)e~™e-^<, as ReC > 0, ImC < 0. 

Set 

MC) = 5(-(C + 1)'/'). (4.145) 

It is now easy to verify that h is analytic and bounded for ReC > 0, and that /i(C) = 
©(e"^!^!) for C — )• ±ioo. By Carlson's theorem, this implies that /i = if a; > 0. Tracing 
back the previous steps, it follows that 5 = 0, M = 0, and = 0, which proves the 
vanishing lemma. □ 
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Remark 4.9 The proof of the vanishing lemma does not apply if either a is not real or 
/5 is not purely imaginary. The first failure is that the function H would not be analytic 
across (0, ^) in this case. A further problem in the proof would be that the matrices 
V\ + and V2 + lose their symmetry, which results in non-zero off-diagonal entries 
in (I4.14ip . It is of course possible that the vanishing lemma can be proven in a different 
way. Another possibility is that, given a and /3, the RH problem is not solvable for 
certain isolated values of x. 

5 Asymptotics for Toeplitz determinants 

Using the identities of Proposition 14.41 and the Fourier representation for V{^z\ we can 
rewrite (|3.46p in the form, with x = 2nt, 



at 



-t 



{a + f3)n - {a^ - /3^) 



+ -a(x) — vix) < a + a [ ; — 

i ^ ^ ^ ' \ \t sinht 



sinht 



kt 



k=l 



k=l 



+ 2 ^ k{Vk + F-fc) cosh(A:t) \ + 0(l/n)$( 



k=l 



(5.1) 



The expressions (fi^Mj) . (fi^HTl) . (|4.118p yield the x expansion for a in (fOB]) : 
and the expressions (I4.22|) . (I4.119P imply the x — )• +00 expansion in ()1.26p . 

Because of the uniformity property of the error term in (jS.ip , the integration of this 
identity from e > to some t < to gives uniformly for any < e < t, 



In Dn{t) = In Dn{e) + {a + P)n{t - e) + 



k=l 



-tk 



Vk-{a + f3)- 



-tk 



00 



kVkV-k + {a-(3)Y, Vke-''' + (a + /?) ^ V-ke 



■ke 



k=l 



+ 



2nt 



aix 



k=l 
,2 n2 



k=l 



dx + (a^ -/3^)ln{n(l 



_J 2ne 



(0^-/32) Inn + /?„(t) + 0(l/n); 



Rn{t) = - v{2nt) jo-ha Q - 4^^^ + 2 ^ /c(Vfc cosh(/!;t)| dt. 



(5.2) 



If a is real and (3 is imaginary, we can take as a path of integration the interval 
[e, t] of the real line as, according to Section 14. 4| the functions (t{x) and v{x) are real 
analytic for positive x. The estimates (jl.26p . (I1.46P ensure integrability at j; = and 
X = +00. In particular, the term in the square brackets in (15. 2p converges if e — ?■ 0. 

For arbitrary /3, Re a > —1/2, we can choose a path of integration and the end- 
point t to avoid possible singular points {xi, . . . , x^}. The estimates (jl.26p . (I1.46P were 
obtained above for positive x. The restriction to real x was only imposed for simplicity 
of notation. In fact, it is easy to verify that the estimates (ll.26p . (I1.46P hold for any 
path to zero and infinity within a sector —tt/2 + 5 < axgx < 7r/2 + 5, Q < 5 < 7r/2. 

We have for Rn{t) in (fO]) : 



\Rn{t)\ <C [ \v{2nu)\du = 0{l/n), 
Jo 



n 



00, < t < to- 



(5.3) 
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Now recall that (|5.2|) is uniform in e, and lnD„(t) is continuous at i = 0. Therefore, 
taking the limit e — >• in ()5.2p and using the Fisher-Hartwig asymptotics p.Sp for 
In Dn (0) gives the expression p.24p of Theorem 11.11 

This concludes the proof of both Theorem 11.11 and Theorem 11.41 
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